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Abstract 

A method for approximating dark soliton solutions of the nonlinear Schrodinger equation under 
the influence of perturbations is presented. The problem is broken into an inner region, where 
core of the soliton resides, and an outer region, which evolves independently of the soliton. It 
is shown that a shelf develops around the soliton which propagates with speed determined by 
the background intensity. Integral relations obtained from the conservation laws of the nonlinear 
Schrodinger equation are used to approximate the shape of the shelf. The analysis is developed for 
both constant and slowly evolving backgrounds. A number of problems are investigated including 
linear and nonlinear damping type perturbations. 

Perturbation theory as applied to solitons which decay at infinity, i.e. so-called bright solitons, has 
been developed over many years cf. [U [2j [3]. The analytical work employs a diverse set of meth- 
ods including perturbations of the inverse scattering transform (1ST), direct multi-scale perturbation 
analysis, perturbations of conserved quantities etc; the analysis applies to a wide range of physi- 
cal problems. In optics, a central equation which describes the envelope of a quasi-monochromatic 
wave-train is the nonlinear Schrodinger (NLS) equation,, which in normalized form is given by 

iU z + ^U tt + n\U\ 2 U = 

where D,n are constant. In this paper we will consider the NLS equation in a typical nonlinear 
optics context where D corresponds to the group-velocity dispersion (GVD), n > is related to the 
the nonlinear index of refraction, z is the direction of propagation and t corresponds to the retarded 
time. In this form, the sign of D determines whether the light focuses or defocuses. In the anomalous 
GVD (or self- focusing nonlinearity) regime the NLS eq. exhibits so-called "bright" solitons which are 
pulses which decay rapidly at infinity. In this case the solitons are formed due to a balance between 
dispersion and self-focusing cubic nonlinearity. 

On the other hand, in the normal GVD (or self-defocusing nonlinearity) regime decaying pulses 
broaden and bright solitons of the NLS eq. do not exist. Instead solitons can be found as localized 
dips in intensity which decay off of a continuous-wave (cw) background. These dark solitons, which 
are termed black when the intensity of the dip goes to zero and grey otherwise, are also associated 
with a rapid change in phase across the pulse. The experimental observations of dark solitons in both 
fiber optics [I] and planar waveguides [5] sparked significant interest in the asymptotic analysis of 
their propagation dating back two decades. 
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Using the above mentioned analytical methods, the propagation of bright solitons under perturbation 
is described by the adiabatic evolution of the soliton parameters; i.e. the soliton's height, velocity, 
position shift and phase shift. However, the non-vanishing boundary of dark solitons introduces serious 
complications when applying the perturbative methods developed for bright solitons. In early work, the 
particular case of added linear loss was studied both numerically [6] and analytically [7] . The analysis 
was specifically for black solitons and solved explicitly for higher order correction terms. These results 
were re-derived [8] with a more straightforward method. The method was extended to grey solitons and 
general perturbations but only for the two of the four main soliton parameters; background height and 
soliton depth, were determined. The evolution of the background was shown to be independent of the 
soliton by Kivshar [9] where the asymptotic behavior at infinity was used to separate the propagation of 
the background magnitude from the rest of the soliton. The amplitude/width of the soliton 'core' was 
determined via a perturbed Hamiltonian. Of the methods proposed many have employed perturbation 
theory based on 1ST theory. In [lOj orthogonality conditions are derived from a complete set of squared 
Jost functions (eigenfunctions of the linearized NLS operator) [11]; from these conditions the soliton 
paramters are, in principle, determined; over the years various corrections/modifications have been 
made to the details [1211151113] . 

In this paper we address a central issue systemic through all these methods. For dark solitons find- 
ing the adiabatic evolution of the soliton parameters (background height, soliton depth, position shift 
and phase shift) alone is insufficient to fully characterize the evolution of a dark soliton. We find 
both analytically and numerically the existence of a shelf which develops around a dark soliton under 
perturbation. The tendency for shelves to generate around dark solitons under external perturbation 
|15] was used to explain discrepancies in the perturbed conservation laws. But, the analytical calcu- 
lation of the core soliton parameters was not obtained. Subsequently the shelf contribution has been 
ignored. However, the shelf is critical in developing the perturbation theory and has a non-trivial 
contribution to the integrals employed to find expressions for the soliton parameters. In this paper, 
we use perturbed conservation laws since they can be easily derived directly from the NLS equation 
and do not require the associated subtleties inherent in the 1ST method and extends to non-integrable 
problems. To carry out the procedure we employ suitable asymptotic information about the higher 
order perturbation terms (beyond the soliton); in general we do not need the exact higher order so- 
lution to solve the leading order problem for the key parameters. Finally we note that shelves in 
soliton perturbation theory have been found earlier in a different class of problems. They were needed 
to effectively understand the KdV equation under perturbation In the KdV eq there is a small shelf 
produced in the wake of the soliton. The height/speed of the soliton, shelf and the additional soliton 
parameter which determines the center of the soliton are all determined by perturbation theory [2J.** 

The outline of this paper is as follows. In Section Q] we pose the problem and illustrate how the back- 
ground evolves under perturbation independent of any localized solitary wave disturbances. Sections 
[2j-[5]set up the basic analysis and a prototypical problem is discussed which helps describe the ideas. 
The method of multiple scales is employed to find the the first order approximation for a black soliton 
under the action of a dissipative perturbation which decays to zero well away from the soliton core. 
The concept of a moving boundary layer is used to bridge the differences between the inner soliton so- 
lution and the outer background. This discrepancy between the approximate soliton solution and the 
background manifests itself as a shelf developing on either side of the soliton. Perturbed conservation 
laws are used to find the growth of the shelf in both magnitude and phase. The analytic results are 
shown to be in agreement with numerical simulations of the perturbed NLS equation. In Sections [6]- [7] 
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the method is extended to grey solitons under general perturbations. Suitable asymptotic information 
about the shelf is obtained from the linear first order perturbation equation; the complete solution 
of the linear problem is not required. In Sections [9]- ?? the perturbation method is applied to some 
physically relevant perturbations: dissipation and two photon absorption. We find that the spatial 
frequency of the soliton differs from that of the background that it resides on. All of the adiabatically 
varying core soliton parameters and the shelf have not been obtained in the many previous studies of 
perturbed dark solitons. In Appendix ?? we derive secularity conditions from Fredholm alternative 
type arguments which agree with the results found from the perturbed conservation laws. 

1 The Boundary at Infinity 

Let us consider the problem of the NLS eq. with normal dispersion D = — 1, n = 1 and with an 
additional small forcing perturbation (we can always rescale NLS to get these unit values) 



where |e| <C 1. Further we will assume a non- vanishing boundary value at infinity; i.e., \U\ -/> as 
t — > ztoo. The effect the perturbation has on the behavior of the solution at infinity is independent 
of any local phenomena such as pulses which do not decay at infinity; i.e. dark solitons. In the case 
of a a continuous wave background, which is relevant to perturbation problems with dark solitons as 
well as in applications to lasers, we have Uu — > as t — > ±oo and the evolution of the background at 
either end U —¥ U (z) is given by the equation 



The above equations completely describe the adiabatic evolution of the background under the influence 
of the perturbation F[u]. Although this is true for all choices of perturbation, we will further restrict 
ourselves to perturbations which maintain the phase symmetry of equation (jl.ip ; i.e., F[U(z, t)e 1 ®] = 
F[U(z, i)]e ld . This is a sufficient condition to keep the magnitude of the background equal on either side 
and a property of most commonly considered perturbations. We assume that at z = 0, u + (0) = u~(0), 
then, since u (z) satisfy the same eq., the evolution is the same for all z. Hence u + {z) = u~{z) = 
u OQ (z). While this restriction is convenient the essentials of the method presented here apply in 
general. The equations for the background evolution (jl.3p can now be further reduced by considering 
the phase difference Acj) QO (z) = (t>toi z ) ~ fiooi 2 ) which is the parameter related to the depth of a dark 
soliton (see below); here ± (z) represents the phase as t — >■ ztoo respectively 




(1.1) 




d 



elm [F[ Uoo \] 



(1.4a) 




'oo 







(1.4b) 
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Thus, while the magnitude of the background evolves adiabatically the phase difference remains 
unaffected by the perturbation. 

Let us now focus on the evolution of a dark soliton under perturbation. To simplify our calculations 
we take out the fast evolution of the background phase 

so equation becomes 

iu z - -u tt + Ol 2 - ul^u = eF[u] (1.6) 
The dark soliton solution to the unperturbed equation is given by 

u s (t, z) = (A + iBt&nh [B (t - Az - to)]) e ia ° (1.7) 

where the core parameters of the soliton: A, B, to, ao are all real, the magnitude of the background is 
{A 2 + B 2 ) 1 / 2 = Uqo and the phase difference across the soliton is 2 tan -1 (^), A ^ 0. When ^4 = 
equation (|1 .T|) describes a black soliton, which has a phase difference of ir. 

Below, we employ the method of multiple scales by introducing a slow scale variable Z = ez with the 
parameters A, B, to and <to being functions of Z. A perturbation series solution for equation (jl.6p is 
assumed 

u = Uo (Z,z,t) + e Ul (Z,z,t)+0{e 2 ) (1.8) 

The first order approximation uo(Z,z,t) should satisfy the slowly varying boundaries from equations 
(|1.4p . which means two of the parameters are already pinned down A[Z) = u OQ (Z) cos ^ 2 °° and 

B(Z) = Uoo(Z) sin ^2°° where we identify o"o = ^ "t^ .. This, however, means that the shape of the 
soliton at first order is determined by non-local effects. This, however, means that the shape of the 
soliton at first order is determined by non-local effects and this suggests that higher order terms might 
be required in order to characterize how the soliton evolves under perturbation. 

2 The First Order Correction 

We write the solution in terms of the amplitude and phase: u = qe 1 ^ where q and are both real 
functions of z and t so equation (jl.6p becomes 

iq z -4> z q-^ (qtt + i^Mt + q(i<t>tt - 4>t)) + (\q\ 2 ~ U L)Q = eF i u ] 

Once we introduce the additional multiple scale variable Z = ez; the real and imaginary parts of the 
above eq. are: 

q z = -(2</Hqt + q<M+e(hn[F[u]]-qz) 
<f> Z q = ~(qtt-(f>tq) + (\<l\ 2 -uloh + eiReiFiu}} - 4> z q) 
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Expanding q and 4> as series in e: q = qo + eq\ + 0(e ) and (j) = <f>o + e4>\ + 0(e ), we have at O(l) 

<7o 2 = ^(2<t>otQat + qo<t>ott) (2.2a) 
</>0z<?0 = —zilott ~ (PotQo) + (\qo\ 2 ~ «L)<?o (2.2b) 
with the general dark soliton solution 



q = (A{Zf + B(Z) 2 t a nh 2 (x)) 1/2 (2.3a) 

+ <J (Z) (2.3b) 



tan 1 



A(Z) 

x = Bit- A(es)ds - t {Z) \ (2.3c) 



o 



For a black soliton the form of solution (12.31) is taken to be 



q (Z,z,t) = u^tanh [u^ (t - t (Z))} (2.4a) 
MZ,z,t) = a (Z) (2.4b) 

where we note that in this representation qo is allowed to be negative. 
At 0(e) we have 

Qiz = ^ [ 2 (^ot<?it + qot<i>it) + qo&ut + qi<t>m\ + Im[F] - q z (2.5a) 



h z qo = -q\\ 



^z- \ [qm ~ (2^Wit)<?o - 4>ltqi] + 3< 7o<7i ~ u L<7i + - (/> z<7o (2.5b) 



Solving for qx,(j>i may be difficult and is not ideal for a general method seeking to obtain the basic 
core parameters of the soliton to require finding its solution. Nevertheless, for clarity we will begin 
with an example in which higher order terms can be found explicitly. This will also provide useful 
insight for the general problem. Consider 

F[u) = iju tt , 7 > (2.6) 

and for concreteness, at a leading order we assume a black pulse which then yields 

qoz = -tozqot (2.7a) 
(poz = ctqz (2-7b) 



If we look for a stationary solution, q\ z = <\>\ z = 0, and note that (pot = (pott = then equations (|2.5p 
reduce to 

= i [2{q 0t (f) lt ) + qo<put] + Im[F] + t oz qot (2.8a) 
= --qm + 3<?o<7i ~ u L<7i + Re[F] - a oz qo (2.8b) 
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where Im[F] = "fqott and Re[F] = 0. 
First we look at equation H2.8aj) 



qot<j>it + -^Qo4>m = - jqou - tozqot 



(2.9) 



which after multiplying by qo, using properties of the leading order solution and integrating, yields 



ht = g 7<7o ~ toz + ciq 



(2.10) 



Since q has a singularity at to we set c\ = to have a bounded solution. After integrating we are 
left with 

4 

01 = -7ln[cosh(uoo(t - t ))] - hzt + c 2 (2.11) 
Asymptotically for large t we have 



Ht = gT^oo - t z, 



(2.12) 



where the superscript ^ indicates the value of a function as t — > ±00 respectively. 



We can also solve explicitly for q\. After a change of variables x = u^t — to) and substituting in for 
qo equation (|2.8bp becomes 



qixx + (6tanh 2 (:r) — 4)qi = — 2— — tanh(a;) 

U 00 

The homogenous problem is now a special form of 

Qxx + (n(n + l)sech 2 x — n 2 )Q = 

which has the bounded solution 

Q = sech n x 

With this solution we can use reduction of order to solve equation (|2.13l) 



qi 



c\ + C2 ( - sinh(4x) + 2sinh(2x) + 3x ) H — ( x — — sinh(4x N 

■ 1 J 8ti„o V 4 



sech 2 (x) 



By looking at the asymptotic behavior at x — > ±00 



1 / &0Z\ ±2x 

qi ~ c 2 e* 

16 V ^00 / 



(2.13) 

(2.14) 
(2.15) 

(2.16) 
(2.17) 



we see that to avoid blow up we must take c% 



. Furthermore, we require that the full solution 



qi vanish at t = to so that the u remains anti-symmetric. Now the unique solution to (|2.8bp is 



[sinh(2uoo(i - to)) + 2uoo(t - t )] sech 2 (u OC) (t - t )) 



Looking at the asymptotic behavior as t — > ±00 we have 



qi -> (2*4,) 1 a oz qo 



&oz_ 

2tioO 



(2.18) 



(2.19) 
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3 Boundary Layer 



Notice that q\ ft and <j)\ ft as t — > ±oo. As a result, the solution u m (qo + eqi)e l ^ 0+e ^ to 
order e does not match the boundary conditions at infinity. Thus, our problem is now broken into 
two areas: the region which matches behavior at infinity and is unaffected by the soliton, and the 
region in which the O(e) correction term is correct. We introduce a boundary layer in which there 
is a transition from a nonzero value in the perturbation term to zero (see also ?]). Note in this 
section we will consider the more general case when Uoo is a function of Z = ez. We find the behavior 
of this boundary layer, where the regions are matched, by looking for a near constant wave solution. 
For this we return to equation (|1.6|) and seek a solution perturbed around the solution at infinity, say 
u ~ (uqo + ew)e 1 ^ +e ®^ where w and 9 are real functions of z and t; the equation is satisfied at O(l) 
and we have at O(e) 



2 Uoo + iw. 



lUr 



At + VHt] + 2n^ 



F[uoo + ew] 



du r 



dZ 



dZ 



(3.1) 



After substituting in equations (II. 3|) and (II. 4h the right hand side is F[uoo + ew] — -Ffuoo] eF[w\. 
Thus, the right hand side is actually a higher order term and may be dropped and as a corollary the 
boundary layer is independent of perturbation. We now break (I3.ip into real and imaginary parts 



w. 



- Uoo V tt 



-w tt 



(3.2a) 
(3.2b) 



Taking a derivative with respect to z of equation ()3.2bp and then substituting in for 6 Z we get 

i r i , i 



It'; 



-wtt tt 



ulo W tt 



1 



rWtttt 



Similarly, taking a derivative with respect to z of equation (|3.2ap and then substituting in for w z we 
get 

1 



Ur, 



Leaving us with 



2ul c w z - -w ztt 



-{-Uoo9tt)tt 



u lo^tt - ~^tttt 



II':: 



ulo w tt ~ jWtttt 



UoaOtt — jOtttt 



(3.3a) 
(3.3b) 
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which is the same equation for both functions, though we will need a different solution for each. This 
is because of the differing boundary conditions to correctly match the inner region to the outer region. 
To the left of the soliton these boundary conditions are 



w(— oo) = w(oo) = g 1 (3.4a) 

0(-oo) = t (oo) = 0i t (3.4b) 

To the right of the soliton these boundary conditions are 

w(— oo) = qf w(oo) = (3.5a) 

t (-oo) = (f)f t 6>(oo) = (3.5b) 



If we let w = e l ( kt+ fo ^( z > k ) dz ) ; then the 'dispersion' relation for equation (|3.3ap is found to be 

u 2 = u 2 O0 (Z)k 2 + \k i (3.6) 

For long waves (k <C 1) we have roughly oo(z, k) ~ ±.u 00 {z)k or w = e %k<yt± ^o n °°( ez )). Thus, we see that 
long wave solutions (i.e. |A;| <C 1) move with instantaneous velocity V(z) = ±u oa (z). This is also true 
for equation (|3.3bp . 

With this in mind, we look for solutions to equations (13. 3|) in a moving frame of reference: x = t — Vz 
and Q = z. 

= Ww^ + (ulo - V 2 )w xx - -w xxxx 

#CC = 2V'0£ B + (u 2 ^ — V 2 )6 XX — -Oxxxx 



And, for V = ±u c 



Wtf = 2Vw c>x - ^w xxxx (3.7a) 
^CC = 2V9£ X — -Oxxxx (3.7b) 



We assume that derivatives with respect to x are small; i.e. long waves. There are several ways to 
balance the terms in equations (|3.T|) (see Appendix A), the optimal one being 

% ^ ~ d xxxx < 1 

leaving us with 

= 2Vw c>x - jW xxxx (3.8a) 

= 2V6q x — -.Oxxxx (3.8b) 
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There are now two similarity solutions which we find to satisfy the boundary conditions (|3.4p and 
(|3.5p derived from matching the two regions. First, by making the transformation 9 X = f and 
£ = x/C 1 / 3 in equation ()3.8bp we get 

= -Vlf - -f" 

which can then be further reduced by the transformation f = f' and £ = —2 (y) 1 ^ 3 £ to get 

= /" - £/ (3.9) 

Equation (|3.9j) is the well known Airy equationwith general solution /(£) = ciAi(£) + C2Bi(£) where 
Ai(£) and Bi(£) are special functions defined in terms of infinite series or improper integrals. Since, 
Bi(£) grows exponentially we take c 2 to be 0. For V = — «oo we are looking for a solution 6 which 
goes to zero as x — >■ — oo and as a direct result / = 9 X — > as £ = x/C}^ — > — oo. With this we can 
now unwrap the transformations made earlier. If we consider the boundary conditions on the left of 
the soliton 

rx [-ax/C, 1 ^ 



6((,x) = c 1 / Ai(s)dsx (3.10a) 



OO J — OO 



1 /3 

where a = —2 (y) and c\ = cf>Q t . Note that the sign of a depends on the sign on V. In the same 
way on the right of the soliton we find that when V = the solution is 

0((,x) = c 2 / Ai(s)dsdx (3.10b) 

J oo J —oo 

This solution matches the boundary conditions for the phase 6 with C4 = (^i. 
To get the other solution we begin by factoring out a derivative with respect to x in equation (|3.8a|) 

= (^Vw^ - ^w xxx ^j 
c 3 = 2Vw{ - -w xxx 

For this to satisfy the zero boundary condition (on either side) it must be C3 = leaving us with 

= 2VW{ - jWxxx 

Which under the same procedure used above has solution 

w(x) = c 4 / Ai(s)ds (3.11) 

J —00 

for both V = —Uoo and V = Uoq. An important point is that there are two boundary layers moving 
away from the soliton solution with speed Uoo generating a shelf. This shelf has a non-negligible 
contribution in the integrals employed in soliton perturbation be it in the form of secularity conditions 
or the conservation laws which we will be employing. 
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4 Perturbed Conservation Laws 



We still need to solve for the slowly evolving parameters cro(Z) and to(Z) for the black soliton. This 
will be done by deriving equations for the growth of the shelf from the perturbed conservations laws 
associated with perturbed NLS (11.6p . The shelf is described by the asymptotic parameters qf 1 and 
4>f t , which are in turn expressed in terms of oqz and t^z- In general we will use the Hamiltonian H, 
the energy E, the momentum /, and the center of energy R. 

2 



H 

E 
I 
R 



—oo 
oo 



1 


du 


2 


dt 



1 / 2 



2\2 



dt 



/oo 
- \U\ 2 ] dt 
-co 

Im [uu[] dt 
t (t& - \u\ 2 ) dt 



oo 
oo 



-oo 
oo 



(4.1a) 

(4.1b) 
(4.1c) 
(4.1d) 



Note that since the standard total energy (Ejr tai = f \u\ 2 dt) would be infinite, we define the energy of a 
dark pulse to be the difference of the total energy and the energy of a continuous wave of corresponding 
magnitude. For unperturbed NLS the first three integrals are conserved quantities while the last can 
be written in term of the momentum, i.e.; ^ = —I. Evolution equations for these integrals may be 
easily obtained from equations (jl.ip and (jl.4p 



— = 2elm 



dH 

dz 
dE 
dz 
dJ_ 
dz 
dR 
dz 



E^-v? + 2eRe 

dz 



F[u]u* z dt 
- F[u]u*dt 



— oo 
oo 



2eRe 
-I + 2elm 



F[u]u* t dt 



t (-Ffuoojiico - F[u]u*) dt 



For our example problem only the energy and the momentum equations will be used. 

5 Full Black Solution 



(4.2a) 
(4.2b) 
(4.2c) 
(4.2d) 



We begin with the perturbed conservation of energy 



d_ 

dz 



/oo 
-oo 



|u| 2 l dt = 2elm / Flu^u 



F\u]u*dt 



Substituting in u = qe 1 ^ ', F[u] = i^/uu, T = t — to and taking the terms up to 0(e) we have 

d ro ° 
dz 



/oo roo 
[qo - u lo + rtqoqi] dT = 2e jq OT TqodT 
-oo J — oo 



(5.1) 



(5.2a) 
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At O(l) equation (jOa) is satisfied: £ f [q% - u 2 ^] dT = 0; At 0(e) we have 

^ j-UaoZ r-oo 

— / q oqi dT = -7 / q 2 0T dT (5.2b) 

J —UooZ J —CO 

This is an equation for the change in energy caused by the propagation of the shelf. Notice that on 
the left hand side of equation (|5.2b|) we are only integrating over T G [— UooZ, Uqqz], the inner region 
around the soliton defined by the boundary layers found in the last section. Since qo and q\ are only 
functions of T, we can apply the fundamental theorem of calculus to arrive at 

., 4 

Uoo [qi(u oo z)q (u oo z) + qi(-u oo z)q (-u oo z)} = "T^oog (5.2c) 

And, for large z (although in practice UqqZ only needs to be modestly larger than the the full-width- 
half-max), we take go ±^oo and qi — > qf leaving us with 

4 

qf - Qi = -^Uool (5. 2d) 

By substituting in the asymptotic approximation (|2.19p found early for qf, we arrive at an expression 
for <7o 

croz = -y^ u lo ( 5 - 3 ) 

Next, we consider the modified conservation of momentum 

d f°° roc 

— Im / uu* t dt = 2eRe / F[u]u*dt (5.4) 

J —oo J —oo 

Again, we let u = qe l< ^, F[u] = i^/uu, T = t — to and use the perturbation expansion for u up to 0(e) 
so that equation (I5.4|) becomes 



d 

dz 



/oo poo 
[<PoTqo + e(2(j) 0T q q 1 + 4>uql)] dT = e2Re / i^qoTTqardT (5.5a) 
-oo J —oo 



which in turn reduces in the same way as the conservation of energy to 

</>+ + 4> u = (5.5b) 

By substituting in the asymptotic approximations (|2.12j) found early for (j)f t , we arrive at an expression 
for to 

toz = (5.6) 

This can now be compared with numerics. The magnitude and phase are depicted in Figure [T] and 
Figure [2] respectively. Here we see the inner region, discussed earlier, is t € (—30, 30) where the 
asymptotic solution matches the numerics and rest of the domain constitutes the outer region where 
the asymptotic solution diverges. The boundary layer shown in Figure [3] compares the solutions 
(13. IIP and (]3.10ap to numerics and illustrates how the inner and outer solutions are connected. The 
propagation of this boundary layer can be seen in Figure S] where contour plot contains the soliton 
extending down the middle and shelf extending out from it. The speed of the boundary layer matches 
the speed predicted by our long wave approximation in section [3l 
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Figure 1: Numerical results plotted against the asymptotic approximation for the magnitude \u\ up 
to O(e). Here z = 30 and e 7 = 0.05. 
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Figure 2: Numerical results plotted against the asymptotic approximation for the phase 4> up to O(e). 
Here z = 30and £7 = 0.05. 
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Figure 3: Asymptotic approximation \u\ ~ Uoo + eto and 
compared to numerics. Here z = 30and 57 = 0.05. 
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Figure 4: The predicted location of the boundary layer tsL = ±^oo^ displayed over a contour plot of 
\u\. Here €7 = 0.05. 
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6 The Grey Soliton 



Consider a grey soliton with velocity A(Z) and A 2 + B 2 = u 2 ^. Let u = qe % ^ where q > and <fi are 
real functions of z and t and introduce moving frame of reference T = t — J Q Z A(es)ds — to and £ = z, 
so equation (jl.6p becomes 

zu^ — zAut — -utt + (M 2 — u to) u = (6-1) 

And, then using u = ge*^ 

KQC + *'0f?) - iA(qT + - - [g TT + i2(f) T q T + (i<p TT - 4>t)o\ + Q 3 ~ u^q = eF[q, 0] (6.2) 

This is now broken into real and imaginary parts 

= Aq T + -(2(f) T q T + #tt) + elm [F [it]] (6.3a) 
( q = A4> T q-^{q TT -cj>lq) + {\q\ 2 -u 2 ^)q + eR%[F[u^ (6.3b) 

We now write equations (|6.3b|) in terms of the slow evolution variable £ = eZ and series expansions 
Q = Qo + e Qi + 0(e 2 ) and <p = 00 + £01 + 0(e 2 ). At 0(1) the equations are satisfied by the soliton 
solution (Gil. 



At 0(e) we have 

qi{ = Aqi T + - [2(4> 0T qiT + (PitQot) + 0itt<?o + 0ott<?i] + Im [F[u ]] - qoz 

4>KQo = -<PocQi + A((j)o T q 1 + far) - ^(<?itt - 0oW - 20 O t<?o0it) + MoVi ~ u^qi 
-Re [F[u ]\ - 4>ozqo 

We look for stationary solutions at 0(e) 

= AqiT + - [2(4> 0T qi T + (f>iTqoT) + 0iTT<7o + 0ott<?i] + Im [F[u ]\ - qoz (6.4a) 

= A(cj) Tqi + (piTQo) - ^(qiTT - (Plr^i ~ 2 ^oTqo4>ir) + ^qi - u^qi (6.4b) 



-Re [F[u ]\ - 4> oz q 



Where 



qoz = ~ (AA Z + BB z t&nh 2 (x)) qo 1 + q 0T -t 0Z ) (6.5a) 

<\>oz = (AB z -BA z )tanh(x)qo 2 + <Po T (^-toz^ +voz (6.5b) 

Unlike the black soliton problem, we do not solve this explicitly. We assume a shelf structure similar 
to the one found in our example problem will develop; this is supported by numerical computations. 
Consider equation (|6.4ap in the limit T — > ±oo using qo — > Uoo and u ooZ = /mFfuoJ yields 



= Aqf T + ^f TT (6.6) 
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We assume q\ tends to a constant with respect to t; i.e. qyr — > as t — > ±00. As a result 4>\tt —> 0. 
Then q\ and 4>it both go to constants as t — > ±00 which corresponds to a shelf developing around the 
soliton. Now, from equation (|6.4b|) in the limit T — > ±00 we get 

A<j>f T + 2 UoD qf = -Re [F[ Uoo ]] / Uoo ± ± — 2 ^ + a oz (6.7) 

We define A^o by 

A^o = 2tan -1 ^^ (6.8a) 

the phase change across the core soliton. his is consistent with the soliton parameters A and B being 
expressed in terms of background magnitude, u^, and phase change, A</>o, 

A = Uoo cos(^p) B = Uoo sin(^p) (6.8b) 



Using equation (II. 4p and substituting in (I6.8bj) on the right hand side equation (16.7P becomes 

A<f>f T + 2 Uoo qf = <\>%± + a z (6.9) 

We recognize the right hand side as the derivative of the phase on either edge of the core soliton (See 
Fig [5]). 



7 Grey Conservation Laws 



Next we use the evolution equations (14.2j) to solve for the shelf parameters qf and cf>f t as well as the 
slow evolution variables A, a^z and to. Note that if we find A, then B = (u 2 ^ — B 2 ) 1 / 2 . The edge 
of the shelf still propagates with velocity V(Z) = u oc (Z), however the speed may now vary in z. In 
terms of the moving frame of reference the boundaries of the shelf are 



SUC) 

Sr(C) 



[uoo(es) + A(es)] ds 



[ [uoo(es) - A(es)] ds 
Jo 



(7.1a) 
(7.1b) 



where Sl and Sr give the position in T of the left and right boundaries of the shelf respectively at £. 
Note that A < for all Z, thus the soliton can not over take the shelf. 

We begin with the evolution equation for the Hamiltonian (|4.2ap . 



d_ 



!| ,2 1/ 2 
2KI +2^- 



2\2 



dt = e (u 



00) z 



/OO f'OG 
[ulc - \u\ 2 ] dt + 2eRe / F[u]u* ( dt (7.2) 
-00 J —00 



Substituting in u = (qo + eqi)e l ^ 0+e ^ and changing variables to the moving frame of reference, we 
have up to O(e) 



d_ 

d( 



/oo roo poo 

(Qot + 0or?o) + («« " 4) 2 dT = ^ K>) z / K " Qo] dt - 4eRe / F[u ]Au* 0T dT (7.3) 
-00 J— 00 J —00 
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where both here and later on uq = qoe 1 ^ . The Hamiltonian is unique among the evolution equations 
(14, 2D in that the contribution of the shelf appears only at 0(e 2 ) or higher and may be ignored. We 
now put in the soliton form (j2.3|) to get 

/oo 
F[u ]u* 0T dT (7.4) 
-oo 

Taking a derivative with respect to Z of the equation u 2 ^ = A 2 + B 2 we get 

= 2AA Z + 2BB Z (7.5) 
which can be used consolidate equations (|7.4j) down to 

/oo 
F[u K T cff, (7.6) 
-oo 

The evolution equations for energy (|5,ip and momentum (|5.4j) both remain the same after transform- 
ing to the moving frame of reference 

d 



d( 



/OO /"OO 
[i& - \u\ 2 ] dT = 2elm / [F[ Uoo ] Uoo - F[u]u*} dT (7.7) 
-oo J —oo 



-oo J — oo 

"OO /■oo 



— Im / uv^dT = 2eRe / F[u]u^dT (7.8) 

J —oo J —oo 

The inner region over which qi and 0i are relevant is T G [5x(£), £r(C)], and outside this region 
ql = (j)\T = 0. At O(l) the equations are satisfied and at O(e) we have 

d r s «(o 



d f° R ^> f°° 
B z - — qoQidT = Im / [Flitoduoo - F[u ]uq] dT (7.9a) 

-2{AB) Z - — / [2</> rWi + 0it9q] dT = 2 Re / F[u }u* 0T dT (7.9b) 

Since the integrands on the left hand side are not functions of £, we can apply the fundamental theorem 
of calculus to arrive at 

/oo 
[F^oojuoo - F[« ]«o] ^ (7.10a) 
-oo 

/oo 
F[u }u* 0T dT (7.10b) 
-oo 

We are left now with the evolution of the center of energy 

-77 t( u oo ~ \u\ 2 )dt = -Im / nn t *dt + 2elm / t (F^jii^ - F[u]u*) dt (7.11) 

J —oo J —oo J —oo 

which after transforming to the moving frame of reference is now 

j /■oo /•£ /■oo 

— (T+ / 4 + t )(t& - \u\ 2 )dT = -Im / uu* T dT 

dQ J-oo JO J-oo 

/CO /•( 
(T+ A + t ) {F[ Uoo } Uoo - F[u]u*) dT 
-co Jo 
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After rearranging some terms we have 



d 



Tint 



\u\ 2 )dT 



/oo 
[• 
-oo 



|n oo ~~ \ u \ \dT + Ir& 



uv?pdT 



(7.12a) 
(7.12b) 
(7.12c) 



/oo roc 
[ul, - \u\ 2 ] dT + 2elm / T (Fjuool^oo — F[u]u*) dT (7.12d) 
-oo J — oo 



Line ()7.12a[) yields 



! )dT = -2 [^(uoo - A)g+ + 5 L ( Uoo + A)q~ 



(7.13) 



The terms on line (|7.12bp are the energy equation (|7.7p and cancel out. The terms on line (|7.12c|) are 
calculated up to 0(e) using the previous results (|7.10p 

E(() = 2B-2[S R {Z)q+ -S L (Z)q^}u oo + eE l (Z) + 0(e 2 ) 
/(C) = -2AB - ul [S R (Z)<f>+ - S R (Z)<p u ] + eh(Z) + 0(e 2 ) 

Noting that and Sl are 0(1/ e) in terms of Z. 
When we put everything together in terms of slow evolution variable Z = e£ we get from (|7.12p 

/oo 
T (F[ Uoo \ Uoo - F[u K) dT t + &E X (Z) + eh(Z) 
-oo 

+ [2noo [^^(noo - A)q+ + S L ( UoD + A)gjf ] + 2^,4 [5^ - S^gf] (7.14) 



This breaks into 0(1) terms 



2 [iSjjg^ + S L q x ] + [5r<^ - S L <j> 1T ] = 



(7.15) 



and 0(e) terms which include higher order energy and momentum terms have not been determined 
The six equations (|6.9p . (|7.6p . (|7. 10a|) . (|7. 10b[) and (|7.15p can now be used to solve for the full set of 
seven parameters qf, 4>u( = <^it)> A , ctq, and to. 
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Uoo = Im^Uoo]] (7.16a) 



dZ 



2B ^A = f^e J F[u ]u 0T dT^ (7.16b) 
B z - Im f F[ Uqo } Uoo - F[u }u* dT + Re[FKo]] ) (7.16c) 

J — oo / 



dZ 



qt = 2 toz + A0oz) / - A) (7.16d) 

qi = ^(aoz-Afozj/iuoo + A) (7.16e) 

</>+ r = -2g+ (7.16f) 

0- T = 2 ?1 - (7.16g) 

B z = (u 00 u ooZ -AAz)/B (7.16h) 

A0oz = (2A5 Z - 2B,4 Z ) /<& (7.16i) 



These equation may now be easily solved from top to bottom. We have added equations (|7.16hp and 
(|7.16ip to the list since it is often better to use these formulation for Bz and A<J)qz rather then working 
out B and Acfto explicitly and then taking derivatives. 

The equations found can also be used to show that the soliton and shelf taken together do not break 
the phase component of the boundary condition. By combining equations (17. 6p and (|7.10bj) we arrive 
at 

2{AB) Z + t4 [(uoo - A)^ + (uoo + A)fi T ] = 4BA Z (7.17) 
which may be rewritten as 

2AB Z - 2BA Z + [MSr) - MSl)} = (7.18) 

If we define (pi as follows 

A^ = MSr)-MSl) (7.19) 

then €(pi is the phase change across the shelf. Substituting this definition along with (|6.8ap into 
equation ()7.18p we arrive at 

Thus, the total phase change across the inner region remains constant, which agrees with our earlier 
result that A^oo (the phase change from — oo to oo) remains constant for all perturbations. 
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Figure 5: Diagram of a perturbed soliton for both the Magnitude and Phase. Here F[u] 
z = 30, €7 = .05 and Ac6 = 4vr/5. 



19 



8 toz and Asymptotic Behavior 



To find the final parameter to we must find a first order solution. We look for series solution to 
equation (jl.ip of the form u = no + eu\ + 0(e 2 ), and at 0(e) we have 

iuiz + (-^ d t + 2 K| 2 - ul^j ui + ni = F[u ] - iu z (8.1) 

or after changing variables to moving frame of reference T = t — J Q A(es)ds — to, z = £ 

iu K + (-iAdr - + 2 K I 2 - u'Cj ui + uj = F[«o] - *«02 (8.2) 

n z = ^ze 4CT + ^ (n - Ae i(T ) + n 0T (-toz + + ia z u (8.3) 



Here 



If we look for stationary solutions, this can be written as a system of coupled second order differential 
equations 

LUi = G[n ] (8.4a) 



where 
and 

L 



Ui=( ft 1 !) G[n ]=f Re |^oj--oz] j (g4b) 

1 Im[mJ / V Im [F[u \ - iu oz \ ' 



±df + (3 A 2 + B 2 tanh(BT) - n^) Ad T + 2ABtanh(BT) 

-Ad T + 2A J Btanh(.BT) - \d\ + (A 2 + S^tanhfsT) - n^) 



3.4c) 



This system has homogeneous solutions 



U 12 = ( Btan f T ) (8.5b) 



U - = < S eJ(BT) ) ^ 
lbJ 

s(srtanh(sr)-i) 

Ui -> 1 A(-BT+§imech 2 (i?T) + §tanh(5T)) 1 ' ' 



3BTsech 2 (BT) + 4tanh( J BT) + tanh(£T) cosh(2£T) 

and using the method of undetermined coefficients we can obtain a particular solution. 

To put u\ in terms of our magnitude and phase functions qo, q\, (j)Q and (fii, we expand our previous 
approximation for u 

u = (qo + egi)e^ 0+£< ^ = q^ +e( gi + ifaqo) e**° + 0(e 2 ) (8.6) 
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so that 



uo = qoe^° (8.7) 

ui = {qi +i<p iqo )e^° (8.8) 

= [qi cos((/) ) - 4>iqo sin(</> )] + i [<?i sin(^ ) + (piQo cos(^ )] (8.9) 

Finally, toz is determined be taking the asymptotic behavior of the solution u% as t — > ±oo 

v£ t = -tf t (±B) + i<ft t (A) (8.10) 



9 Example for Grey Solitons 

Let us return to the perturbation F[u] = i'juu, however we now consider the evolution of a general 
dark soliton with Uoo(0) = 1. As was the case for black solitons, the background height is found 
to be constant from equations (II. 4p . In Fig [6] we see that the velocity of the soliton does not effect the 
velocity of shelf which still moves with velocity V = ±Uoo- Using the equations derived in sections [6] 
and [7] we now solve for all relevant parameters 



A z = 





(9.1a) 


&0Z = 


4 2 ■ 3 / n 

-g7«ooSin (a) 


(9.1b) 


qt = 


2 

-g7(«oo + A) sin (a) 


(9.1c) 


q\ = 


2 

-g7(«oo - A) sin (a) 


(9.1d) 


ti* = 


4 

^7 («oo + -4) sin (a) 


(9.1e) 


<t>u = 


4 

-g7C«°o -^4) sin (a) 


(9.1f) 



where a = and Ac^o is the phase change across the core soliton as defined in equation (16. 8aj) . 

Unlike the speed of the shelf, the magnitude of the shelf does depend on the soliton's velocity (which 
is in turn related to the soliton's depth and the phase across the soliton). As illustrated in Fig [7] the 
shelf grows shallower behind the soliton for lager speeds (or smaller phase change A0o). The extra 
phase <jq(z) = — ez^ju^ sin 3 (a) induced by the perturbation means that the spatial frequency of 
the soliton is different then the frequency of the cw background that it lies on. Though do evolves 
adiabatically the soliton eventually becomes noticeably out of phase from the background as shown in 
Fig [8] Here the background phase (<p + and is constant since the fast evolution of the background 
phase was take out in equation (jl.ip . 

10 Conclusion 

In conclusion, we have derived a new approach to dark soliton perturbation theory which break the 
problem into an inner region around the soliton and an outer region equal to the boundary at infinity. 
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Figure 6: The predicted shelf edge overlaid on the contour plot of numerical results. Here 67 = .05 , 
and A0o = 47r/5. 
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Figure 7: The shelf height found numerical for various values of A0o plotted over the asymptotic 
approximations. Here Uqo = 1 and = .05. 
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Figure 8: ctq(Z) plotted against the phase at plus and minus infinity along with the phase at the center 
of the soliton. Here Mqo = 1, 67 = .05 , and A</>o = 47r/5. 



We find that under perturbation a dark soliton develops a shelf the edge of which propagates out 
at speed equal to the magnitude of the cw background. Analytically this shelf arises due to the 
difference between the perturbed soliton and the soliton solution of unperturbed NLS which satisfies 
the boundary conditions. It was shown to be possible for the soliton to have a different frequency 
then the cw background. The method extends to general perturbations and works for both moving 
and constant background. For example perturbations the asymptotic approximation calculated was 
compared to numerical results. These comparisons confirm the existence of the analytically predicted 
shelf and supported our claim that the adiabatic approach alone is insufficient to fully describe the 
dark soliton evolution. The main conclusion is that the non-vanishing background and soliton must 
be treated separately to have a consistent perturbation theory for dark solitons. 
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